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I. INTRODUCTION 


: This final report gives details of the second year of a 
four year numerical ocean circulation modeling program for the 
Gulf of Mexico. The aim of the program is to progressively 
upgrade in modest increments an existing numerical ocean circu- 
lation model of the Gulf so that the final model has a hori- 
zontal resolution of about 10 km and vertical resolution 
approaching 1 to 10 m in the mixed layer, 10m at the thermo- 
cline and 100 m in the deep water. Throughout the four year 
period, the validity of the upgraded model will be continuously 
tested, and velocity field time series delivered periodically 
based on the most realistic simulation of Gulf circulation 
available. 


Experiments in the first year were with the existing 
NORDA/JAYCOR two layer hydrodynamic primitive equation ocean 
circulation model of the Gulf on a 0.2 degree grid. They 
concentrated on correctly specifying the coastline and bottom 
topography for maximum realism in circulation simulation, and 
on how best to includ2 wind forcing. 


Experiments in the second year were also with the 
existing NORDA/JAYCOR hydrodynamic model, but most were on a 
0.1 degree, rather than 0.2 degree, grid. Wind-forced only, 
port-forced only, and wind- plus port-forced simulations were 
generated. Details of selected experiments are presented in 
this report. 


Simulated surface currents sampled every three days for 
more than 10 years were delivered to MMS from Experiment 
201/16.0, which is driven by both wind and port forcing, as 
representing the best simulation available from the second year 
effort. Simulated surface currents sampled every three days 
for two eddy cycles were also delivered to MMS from Experiment 
201/17.0, as representing the best simulation available with 
port forcing only. 


Model parameters for all referenced experiments are 
included in Appendix A. 


II. SUMMARY OF YEAR 1 


The ocean model contains many innovative features, and 
the reader is referred to its detailed discussion in Hurlburt 
and Thompson (1980). In particular, Section 2 (pp. 1613-1614) 
gives the model equations and Appendix B (pp. 1647-1650) 

escribes the numerical design of the model. Since that time, 
the capability to handle general basin geometry has been added 
but this does not affect the description in any major way. . 
Wind forcing is treated identically to interfacial and bottom 
Stress terms, i.e., wind stress appears directly as an additive 
term in the momentum equation (see p. 1614 of Hurlburt and 
Thompson (1980)). 


In terms of ‘realism' Experiment 9 was the most success- 
ful Gulf of Mexico numerical simulation prior to the start of 
this project. The model was driven from rest to statistical 
equilibrium solely by a steady inflow through the Yucatan 
Straits which was compensated by outflow through the Florida 
Straits. Figure 1 compares ‘instantaneous’ upper ocean flow 
patterns just before an eddy is shed by the Loop current (a) 
from the numerical model and (b) from observations by Leipper 
(1970). ‘the ability of the model to simulate observed features 
is clearly demonstrated by this comparison, which is remarkable 
given the simplicity of the model forcing; however, some 
discrepancies remain. For example, the eddy has not penetrated 
as far into the Gulf and is more intense than that shown in the 
observations. Waves can be seen moving around the wall of the 
Loop Current in both the model and the observations, but in the 
model they are at the limit of resolution and therefore 
unrealistically large. Moreover, in the Gulf the waves are 
more pronounced on the eastern wall of the Loop and can form 
strong cold intrusions that may contribute to the eddy shedding 
process (Vukovitch and Maul, 1984). As shed eddies propagate 
westward (Figure 2a) the model spontaneously develops a 
counter-rotating vortex pair (Figure 2b), a structure 
repeatedly observed in the Western Gulf (Figure 3). The roles 
of the wind and the Loop Current eddies in the formation of 
this structure have been a matter of some controversy {Merrell 
and Morrison, 1981). Although wind forcing was not present in 
this simulation, a major role for winds has not been ruled 
out. After spin-up the experiment sheds an eddy once every 390 
days and the observed eddy shedding cycles are very similar. 


There is good observational evidence that the total 
transport Cyrougn , the Yucatan Straits is about 30 Sverdrups 
(Sv) (106 m3 sec-!), but there are very little data 
available on the vertical distribution of the flow. Experiment 
9 allocated 26 Sv to the upper layer and 4 Sv to the lower. 
Allocating less transport to the upper layer would give rise to 
smaller Loop Current eddies, but exactly what range of eddy 
sizes is realistic is hard to quantify. However, there is one 
source of Gulf-wide data that can be used as a guide: maps of 
sea surface variability for the Gulf have been produced from 
all hydrographic, STD and XBT data (Maul and Herman, 1984), and 
from satellite altimeter crossovers (Marsh, et. al., 1984). 

The 20 Sverdrups upper and 10 Sverdrups lower layer distribu- 
tion of inflow transport in Experiment 40 gave rise to a vari- 
ability ap very similar to that obtained from the satellite 
(Figure 4); these maps agree more closely with each other than 
with the map from hydrographic data (Figure 5). Based on the 
agreement of variability maps, the mean sea surface from Exper- 
iment 40 (Figure 6) may well be the best mean available for the 
Gulf (Thompson, 1986). Experiment 40 became the baseline 
port-forced experiment from the Year 1 project. The model 
parameters were: 


e upper layer inflow transport = 20 x 106 m3 
sec-l (= 20 Sv), 


@ lower layer inflow transport = 10 x 106 m3 
sec-l (= 10 Sv), 


@ wind stress = 0, 

@ horizontal eddy viscosity, A = 300 m2 sec! 

e grid spacing = 20 by 22 km (0.2 by 0.2 degrees), 
@ upper layer reference thickness, Hl = 200m 

@ lower layer reference thickness, H2 = 3300m, 

e minimum depth of bottom topography = 500m, 

@ beta, df/dy = 2 x 1071! m-! sec-l, 


e Coriolis parameter at the southern boundary, £ = 4.5 
x 1075 sec -1, 


e Gravitational acceleration, g = 9.8m sec~2, 
e reduced gravity, g' = .03(H1+H2)/H2 m sec~2, 
® interfacial stress = 0, 


® coefficient of quadratic bottom stress = .003; 


e time step = 1.5 hours. 


Experiment 68 is identical to experiment 40 except for 
the addition of wind-forcing after the port-forced circulation 
has fully spun up. The winds used were from the Navy Corrected 
Geostrophic Wind data set for the Gulf of Mexico which is based 
on the Navy's twelve hourly surface pressure analysis from 1967 
to 1982. Wind input was every 12 hours; at first 1967 winds 
were used repeatedly to spin-up the wind-driven flow and then 
winds from 1967 to 1977 were applied in sequence. Figures 7 to 
10 show upper layer currents (i.e., vertically averaged 
currents above the thermocline) every 60 to 120 days for 300 
days. Vectors are only drawn at every second point (i.e., 
every 0.4 decrees) to improve readability. Figure 7 shows the 
furthest norchward penetration of the Loop Current ever 
attained by the ocean model, similar configurations have been 
seen in the Gulf. After the eddy breaks off, the Loop Current 
intrudes onto the Florida Shelf and some of the flow splits off 
to the north for a brief time (Figure 8). Similar intrusions 
have been observed in the Gulf, but the model's inadequate 
representation of shelf topography makes it likely that the 
simulated currents in shallow areas (say less than 100 m) are 
too high. A persistent anticyclonic gyre in the northwest Gulf 
has been a feature of almost all Gulf simulations performed to 
date. The addition of wind forcing in Experiment 68 has 
increased its average size and its effect on incoming Loop 
Current eddies (Figure 9). The presence of a gyre in this 
position is explainable by the northward migration of anti- 
cyclonic eddies along the coast of Mexico until the continental 
Slope bends eastward and they can go no further; the effect is 
magnified because the winds also tend to produce an anticyclon- 
ic gyre at the same spot. However, in the Gulf the gyre 
probably dissipates relatively rapidly against the shallow 
shelf area. 


Simulated surface currents from Experiment 68 have been 
delivered to MMS as representing the best simulation date 
available at the end of Year 1. They consist of velocity com- 
ponent (u and v) fields on a 0.2 degree rectangular grid 
covering the Gulf area, sampled every three days for 3780 model 
days (10.3 years). 


III. YEAR 2 WIND-FORCED EXPERIMENTS 


Simulations forced by winds based on a seasonal clima- 
tology from ship observations (Elliott, 1979), see Figures 
ll(c,d) to 14(c,d), had been performed before this project 
started, and will be briefly described below. Such wind fields 
are not ideal for driving ocean models since they contain very 
little of the total wind variability and mean wind strengths 
are in general far weaker than instantaneous winds. Recogniz- 
ing this deficiency, NORDA funded JAYCOR to produce a wind set 
for the Gulf based on the Navy's twelve-hourly surface pressure 
analysis, which is available from 1967 to 1982 (Rhodes et. al., 
1986). The geostrophic winds, corrected geostrophic winds, and 
wind stresses (all on a one degree grid covering the Gulf) 
every 12 hours from 1967 to 1982 are on magnetic tape. These 
are available through the MMS Gulf of Mexico Regional Office. 


Figure 15 shows the wind stress and wind stress curl from 
this data set for 0000 and 1200 GMT on 14 January and 0000 GMT 
On 15 January 1976. There is large temporal variability of the 
wind field during this period, as general easterly flow gives 
way to strong northerly flow after a frontal passage in just a 
24 hour period. The wind stress curl field also shows the 
rapid change, from a relatively weak field to a very strong 
field with strong horizontal gradients. Figure 16 shows 
Similar plots for 14 and 15 July 1976. Even in the summer, 
when flow is generally weaker, very significant differences can 
be seen in a short time period. These strong variations and 
very rapid changes in the wind field indicate why the modeling 
of Gulf circulation requires wind data on short temporal scales. 


Figures ll(a,b) to 14(a,b) show the seasonal climatolo- 
gies averaged over the period 1967-1982. The wind stress and 
wind stress curls are much stronger in the winter season than 
the summer season as would be expected. There are persistent 
areas of positive curl over the Yucatan and negative curl in 
the southwest Gulf that are present for all seasons, which have 
not been seen in any previous study of Gulf wind stresses. 
Although not present at all time periods (Figures 15 and 16), 
rey are also the dominant features of the instantaneous curl 

elds. 


Experiment 31 was performed before this project started. 
It is identical to Experiment 40 except that it has a slightly 
smaller coefficient of quadratic bottom stress (0.002 rather 
than 0.003) and is wind forced by the seasonal climetology from 
Ship observations shown in Figures 11 to 14 with no port 
forcing. After spinup the experiment exhibits an annual repeat 
cycle; Figures 17(a) to 20(a) show the interface deviation for 
each season with a 12.5 meter (rather than the more usual 25 
meter) contour interval. 


Experiment 202/11.0 was performed in year 2 of the 
project. It uses a bottom topography that is (before smooth- 
ing) a projection of the 0.1 degree topography onto a 0.2 
degree grid. It all other respects it is similar to experiment 
31 except that the initial wind forcing is the monthly clima- 
tology from the Navy Corrected Geostrophic Wind data set, i.e., 
it uses the wind stresses in Figure ll(c) to 14(c) except there 
are twelve monthly fields instead of fcur seasonal (three- 
monthly) fields. As in all wind-forced cases, the winds are 
linearly interpolated in time between wind inputs so that the 
winds change slightly at each timestep and the model is not 
shocked by a large sudden change in wind stress. After spin up 
this experiment also exhibits an annual repeat cycle, Figures 
17(b) to 20(b) show the interface deviation for each season 
with a 12.5 meter contour interval. The two wind sets 
(seasonal ship winds in 31; monthly modeled winds in 202/11.0) 
give rise to qualitatively similar circulation patterns every- 
where except in the southwest Gulf in the winter and spring. 
However, 202/11.0 has interface deviations that are about two 
to three times greater than experiment 31, which is consistent 
with the generally higher wind stress and wind stress curl 
values in the corrected geostrophic wind data set. In general 
ship observation-based wind fields tend to underestimate wind 
strengths, moreover the geostrophic winds have been corrected 
to agree closely with winds reported from the three National 
Data Buoy Center (NDBC) buoys in the central Gulf. Therefore, 
the circulation strength in 202/11.0 is likely to be more real- 
istic than in experiment 31. In the southwest Gulf experiment 
31 has very strong upwelling in the winter and spring, but 
experiment 220/11.0 exhibits relatively weaker downwelling at 
this time. Again, this is consistent with the very different 
wind stress curl patterns of the two data sets in the region, 
but in this case it is difficult to judge which is more nearly 
correct because independent data are very sparse in this 
region. In the northwest Gulf both experiments show a wind- 
induced west Texas coastal current, as was predicted by Sturges 
and Blaha (1976) from calculations of ‘wind stress curl based on 
Wind data tabulated by Hellerman (1967). Note that Loop 
Current eddy remnants often end up at about the same place as 
the wind-induced anticyclonic gyre, and it would be very diffi- 
cult to distinguish some phases of wind-driven and Loop 
Current-driven circulation in the southwest Gulf without a 
fairly long time history of the feature in question. 


Experiment 202/11.0 was continued after it reached a 
quarci-steady state with monthly iveraged w'nds for each year, 
i.@., 12 wind sets per year, starting in 1967. Figures 21 to 
24 show interface deviation plots for each season of 1971 and 
1972. The effects of interannual wind variability are clearly 
seen, however the strength of the circulation does not appear 
to be significantly different than for the monthly climatology. 


Unlike port-forced experiments, simulations with wind 
forcing only do not provide realistic Gulf wide circulation 
patterns. This is because wind induced flows in the Eastern 
Gulf, except over the continental shelf, are very small rela- 
tive to the dominant, port-forced, Loop Current system. How- 
ever, in tive Western Gulf wind forcing is very significant, 
perhaps as important as the large Loop Current eddies that 
Sigrat< westward from the Loop Current. The simulations with 
bo.’ port ani wind forcing, described later, demonstrate that 
the «ffects of winds and Loop Current are not additive; for 
exampic, the Loop Current variability is significantly 
increased over port-forced only cases, despite the relative 
weakness of winds in the Eastern Gulf. 


IV. YEAR 2 PORT-FORCED EXPERIMENTS 


The bottom topography used in all 0.1 degree simulations 
is shown in Figure 25. It is derived from the SYNBAPS data 
base (Van Wyckhouse, 1973), but there is a flat shelf between 
the 500-m isobath and the coast and the topography has been 
smoothed by two applications of a 9-point real smoother to 
filter out features that could cause 2-grid point numerical 
noise in the model. One advantage of the 0.1 degree grid over 
the 0.2 degree version is that the required smoothing has 
Significantly less effect on the steepness of the continental 
slope in the finer grid model. The model boundary is approxi- 
mately at the position of the original 10-m isobath, but the 
position of the plotted coastline in all figures follows the 
actual coastline as closely as possible, i.e., the model 
boundary is often several points seaward of the plotted coast- 
line boundary. The section of the Caribbean shown is treated 
as land by the model, and the position of the inflow port is 
marked by the termination of ccntour lines in the Yucatan 
Strait. In the 0.1 degree model the inflow port is 16 grid 
points wide, compared to % grid points in the 0.2 degree 
model. This allows the inflow velocity profile to be better 
resolved, i.e., the velocities at the core of the inflow jet 
can be higher in the 0.1 degree case for the same inflow trans- 
port. In all cases described here, the upper layer inflow is 
exactly normal to the port and has a cubic velocity profile 
with the maximum about one third of the way from the west end 
of the port and zero velocities at both ends, and the lower 
layer inflow ‘* exactly normal to the port and has a linear 
velocity profile with the velocity at the western end of the 
port twice that at the eastern end. In last year's 0.2 degree 
experiments, such as experiments 40 and 68, the profile in both 
layers was parabolic with the maximum at the center and zero 
velocites at both ends of the port. In all exyeriments to date 
the inflow transport has been spun up from rest to its required 
value of 30 Sv in the first 300 days and has then been held 
constant for the duration of the simulation. A mean transport 
of about 30 Sv, with a seasonal cycle of about 4 Sv, and with 
highly energetic variations on shorter time scales, is the 
generally accepted characterization of the transport through 
the Florida Straits (Niiler and Richardson, 1973; Brooks, 1979, 
Schott and Zantopp, 1985). A further advantage of the 0.1 
degree grid model is that it allows the use of a lower horizon- 
tal eddy viscosity. Eddy viscosity is used to paramatarize 
Ssub-grid scale processes, i.e., in practice to damp out all 
scales whith are too small to be represented by the model's 
finite difference dynamics. Eddy viscosity can be lower on a 
fine grid because smaller scales are resolved. It is usuai to 
run several simulations with different eddy viscosities (with 
the largest value two or three times the smallest), because 
very low values might allow significant circulation features 
that are only marginally resolved by the model dynamics. 


Experiment 201/13.0 is the baseline port-forced 0.1 
degree experiment. The model parameters are very similar to 
the 0.2 degree, experiment 40, except for the new inflow 
profiles and the lower eddy viscosity made possible by the 
finer grid. The parameters are: 


@® upper layer inflow transport = 20 x 166 m3 
sec™+ (= 20 Sv), 


@ lower layer inflow transport = 10 x 106 m3 
sec~+ (= 10 Sv), 


@ wind stress = 0, 

® horizontal eddy viscosity, A = 100 m2 sec~! 

e grid spacing = 10 by 11 km (0.1 by 0.1 degrees), 

e upper layer reference thickness, Hl = 200m, 

e lower layer reference thickness, H2 = 3450m, 

® minimum depth of bottom topoyraphy = 500m, 

@ beta, df/dy = 2 x 10741 m-) sec-l, . 

e Coriolis paraneter at the southern boundary, f£ = 4.5 
x 10-2 sec +, 

e gravitational acceleration, g = 9.8m sec™¢, 

@ reduced gravity, g' = .03 m sec~2, 


e interfacial stress = 0, 


® coefficient of quadratic bottom stress = .002, 


e time step = 1.0 hour. 


Figures 26 to 28 show upper layer velocities every 90 
days for a complete eddy cycle, as is usual for port-forced 
(only) cases all eddy cycles are similar after spin up. In 
these and all subsequent current plots the velocity is shown 
every 0.2 degrees, i.e., only half the model grid resolution is 
plotted, and qui currents higher than 50 cm sec™+ are drawn 

as 50 cm sec™+, In areas where the current is low and in the 


same direction over several grid points the plotted vectors are 
further thinned by a factor of two, this situation occurs 
infrequently in upper layer plots but is noticeable in the 
lower layer. In Figure 26(a) the Loop Current has repenetrated 
into the Gulf after shedding an eddy about 120 days 

previously. This anticyclonic eddy has reached the coast of 


Mexico and has split into two, and is bracketed by two cyc‘onic 
gyres with the gyre in the southwest Gulf being particularly 
large. There is a cyclonic eddy on the wall of the Loop . 
Current over the Florida Shelf; if the model included the true 
Shelf depth this feature would probably dissipate very 

rapidly. In Figure 26(b), 90 days later, the eddies in the 
western Gulf are smaller, the Loop Current is penetrating 
further into the Gulf and is encroaching less onto the Florida 
Shelf. In Figure 27(a) the Loop Current is about to break off 
an eddy; the previous Loop Current eddy is now quite small and 
has started to move north along the coast of Mexico. In Figure 
27(b) the eddy has detached from the Loop Current, surrounded 
by smaller cyclonic eddies; the previous eddy has migrated to 
the Texas coast and is very small. In Figure 28(a) the Loop 
Current eddy has migrated westward and is about to split into 
two, the Loop Current has started to repenetrate into the Gulf, 
and there is a cyclonic eddy just above it. In Figure 28(b) 
the situation is similar to Figure 26(a) and the cycle is com- 
plete, the small remnant of the previous Loop Current eddy 
appears to have been absorbed into the larger of the two pieces 
of the latest eddy. 


Experiment 201/17.0 is similar to experiment 201/13.0 
except that it has a lower eddy viscosity (A=50 m2 sec~+), 
and a larger coefficient of quadratic bottom friction (0.003). 
Simulated surface currents from this experiment, sampled every 
3 days for 2 complete eddy-shedding cycles, have been delivered 
to MMS as representing the best port-forced experiment to 
date. Figures 29 to 34 show upper and lower layer currents 
every 90 days for one complete eddy cycle. In Figure 29 an 
anticyclonic eddy has just detached from the Loop Current, and 
it has a large associated anticyclonic deep water eddy at its 
leading edge. Figure 30, 90 days later, shows that the Loop 
Current eddy moving into the western Gulf has generated a 
counter-rotating vortex pair in the lower layer, with the anti- 
cyclone leading and the cyclone trailing the upper layer eddy. 
In Figure 31 the Loop Current eddy has reached the Mexican 
coast, and the leading deep water anticyclone has dissipated 
against the continental slope. The deep water cyclone has 
strengthened and a cyclone has formed in the upper layer, i.é., 
an anticyclone in the upper layer has become, in 90 days, a 
cyclone in the lower layer with a vortex pair in the upper 
layer. The Loop Current is repenetrating into the Gulf, and 
there are some signs of a northward counter current in the 
lower layer along the Florida slope. In Figure 32 the old Loop 
Current eddy system (two anticyclonic and one cyclonic eddy in 
the upper layer and a cyclonic eddy in the lower layer) is 
dissipating on the Mexican coast, and the leading anticyclonic 
component has started to move northward. A new Loop Current 
eddy is forming and small cyclonic meanders can be seen moving 
around the wall of the Loop Current. In the lower layer the 
northward counter current has intensified as part of a cyclonic 
Circulation below the Loop Current. In Figure 33 a new Loop 
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Current eddy has just separated from the Loop Current. In 
Figure 34 it is moving into the central Gulf, but the previous 
Loop Current eddy (that formed more than 400 days previously) 
is still strongly in evidence to its southwest. The longevity 
of the first Loop Current eddy appears to be related to its 
breaking into two parts and then recoalescing. 


Port-forced experiments provide realistic Gulf wide 
Circulation patterns. However, there is usually very little 
variability from one Loop Current eddy cycle to the next. The 
Loop Current eddy shedding period is usually approximately 
constant for each experiment, although it varies from experi- 
ment to experiment. Almost all port-forced simulations have 
Loop Current eddies following the same south of west path 
across the Gulf. This is probably the preferred path for 
actual Loop Current eddies in the Gulf. Overall, port-forced 
Simulations provide a very good realization of a typical Loop 
Current eddy cycle, and they demonstrate almost everything that 
was known about Gulf circulation up until the mid to late 
1970's. Since then, however, the observational data base has 
been greatly expanded, and it is now recognized that the Gulf 
circulation has far more variability than these simulations 
allow. For example, the Loop Current eddy shedding period can 
be anything from 6 to 18 months, and Loop Current eddies can 
tg a westerly, rather than south westerly, path across the 
Gulf. 
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V. YEAR 2 PORT AND WIND FORCED EXPERIMENTS 


Great difficulty was experienced in obtaining a long term 
Simulation with both port and wind forcing. The basic strategy 
was as for experiment 68 of Year 1, i.e., spin up to 
statistical equilibrium with port forcing only and then to add 
wind forcing.: As in Year 1 it was found that the port-forced 
Simulation with the lowest eddy viscosity, in this case experi- 
ment 201/17.0, could not accept wind forcing. The reason is 
that the range of possible eddy viscosities depends on the 
maximum current speed, which is increased by about 50% when 
winds are added. Therefor, experiment 201/13.0 was selected as 
the base port-forced experiment. But even in this case the 
addition of wind forcing gave rise to strong upwelling in the 
very small region between an anticyclonic eddy and the coast of 
Texas. This upwelling was so strong that the layer interface 
Surfaced, thus stopping the simulation. The usual method for 
controlling such episodic upwelling in the hydrodynamic model 
is to introduce layer thickness-dependent interfacial friction, 
i.e., friction that varies in time and space depending on the 
instantaneous upper layer thickness. For example, if the 
e-folding time (the time to reduce velocities by half) of the 
friction were set at 16 years for a layer thickness of 200 m, 
then the e-folding tine would be 1 year (still quite small) 
anywhere with a layer thickness of 100 m, but it would decrease 
to 22 days for a 50-m layer, and to only 1.5 days for a 25-m 
thick layer. The basic idea is that the friction should have 
little effect providing the layer thickness stays in its normal 
range, say from 100 m to 400 m, but be sufficiently strong to 
prevent strong upwelling events to surface. A more correct way 
to handle upwelling is to add thermodynamics to the model. 
Upwelling events can then effectively ‘surface’ (without the 
layer interface's surfacing) by entraining mass from the lower 
layer until the densities of the two layers become equal. 
Thermodynamics will be added to the model in Year 3 of this 
program. 


Unfortunately, to perform the subject simulation, experi- 
ment 201/16.0, the interfacial friction's e-folding time had to 
be set to 115 days at 200 m (7 days at 100 m) to prevent 
surfacing, and this is sufficiently strong to modify the 
evolution of desired upwelling events such as the cyclonic 
eddies that move around the wall of the Loop Current. Experi- 
ment 201/16.0 applies monthly averaged winds from 1967, 1968, 
and so on, to experiment 201/13.0 after 6 years of port-forced 
Spin up. It has an interfacial friction e-folding time that is 
115 days at a layer thickness of 200 m in the far west of the 
basin, but this increases to an e-foiding time of about 3 years 
for the same thickness in the eastern Gulf. Figures 35 to 39 
show upper layer currents every 90 days for two Loop Current 
eddy shedding cycles with winds from 1968 and 1969. The 


12 


maximum currents lie between 110 and 165 cm sec~l, against 95 
to 110 cm sec~+ for the port-forced experiment 201/13.0, see 
Figures 26 to 29. The two, ecildy-shedding cycles take about 360 
and 450 days respectively and overall there is more variability 
than in the port only case. The anticyclonic gyre off the 
Texas coast is generally larger than in 201/13.0, as would be 
expected from the presence of a similarly positioned gyre in 
the wind forced case, see Figures 21 to 24. Simulated surface 
currents frou experiment 201/16.0, sampled every 3 days for 
more than 10 years, have been delivered to MMS as representing 
the best Gulf simulation to date. 
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VI. LAYERED MODELS AND BOTTOM TOPOGRAPHY 


Bottom topography is very important for realistic simu- 
lations of Gulf circulation. Although the NORDA/JAYCOR model 
has an excellent representation of bottom topography, it still 
suffers from the fundamental drawback of layered formulations: 
topography must be confined to the lowest layer at all times. 

A layered formulation was chosen for the ocean model because 
alternatives that represent shelf topography more accurately do 
so by increasing the vertical resolution by a factor of 5 to 
10. This leads to corresponding increases in computer time and 
memory that make realistic simulations of Gulf circulation 
impractical (i.e., prohibatively expensive) using such models 
on even the fastest existing supercomputers. 


In all of JAYCOR's two layer Gulf of Mexico simulations 
to date the minimum depth of the bottom topography is set at 
500 m. The layer interface is initially at 200 m, and can 
deepen at the center of Loop Current eddies to about 450 m, so 
there is no danger of the interface intersecting the 
topography. Figure 25 shows the topography used by the 0.1 
degree model. There are broad areas near the coast where the 
model topography has been unrealistically set to 500 m. The 
experiments performed to date demonstrate that topography 
limited to 500 m is sufficient to realistically simulate 
Circulation in the deep water of the Gulf. For example, Loop 
Current eddies show evidence of topographic steering, even 
though there is no direct interaction between the eddy in tae 
upper layer and the topography in the lower layer. 


Continental slope and shelf areas cannot be adequately 
represented by existing layered models since the shelf water 
should be entirely from the upper layer. Apart from the 
Obvious problems this causes for flow on the shelf, the arti- 
ficial deepening of the shelf in these models has at least two 
possibly unrealistic effects on deep water circulation. It 
allows deep water phenomena, such as mesoscale eddies, to 
encroach on ‘shallow' shelf areas, and it drastically reduces 
the potentially very large energy dissipation over the shelf. 
In particular, the Loop Current tends to encroach on the West 
Florida Shelf in the model (although limited encroachments have 
been observed in the Gulf, the present model cannot be expected 
to simulate them realistically), and a Loop Current eddy 
remnant off the Texas coast which is a semi-permanent feature 
of all simulations to date would almost certainly dissipate, or 
move further offshore, in a model that had a better represen- 
tation of the shelf area. 


In order to quantify how well the ocean model is perform- 
ing in shelf areas a comparison has been made, by Evans Waddell 
at Science Applications International Corporation, between 
three years of hydrograhic data obtained at fixed moorings on 
the Florida Shelf as part of a MMS sponsored Gulf of Mexico 
Physical Oceanography Study (Waddell, 1986) and 1080 days of 
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3-YEAR JAYOOR MODEL SIMA ATION 
& OBSERVED DATA STATISTICS 
UPPER LAYER 
DST. WATER PRINCLPAL AXES 
D UT WONG )«=s«éEPTH«s«CDEPTH«U-EAN UA «—sU-MIN U-VARIANCE «V-MEAN «V-MAK «-MIN «= V-VARIANCE «(DEGREES TRUE) 
Max! 25:48.00  84:12.00 100 505 26.6 81.0 -19.3 758.9 -23.9 18.5 -77.9 567.3 130.0 
MiC2 25:53.20 84:19.20 100 180 0.6 36.8 27.8 53.7 4.1 33.0 -61.8 224.6 176.0 
hUX! 25:48.00 84 : 12,00 100 505 26.6 81 0 ~19.3 758.9 23.9 18.5 -77.9 567.3 130.0 
tWE2 27:25.20 684:37.50 100 180 0.5 20.8 -26.5 18.7 3.2 %.1 43.5 104.8 179.0 
LA 26:12.00  83:00,00 100 500 3.2 37.7 -M4.4 139.3 3.5 47.2 -22.0 118.4 129.0 
MYF2 26:14.50  83:13.30 40 50 0.0 13.7 29.8 6.0 -1.0 21.1 20.1 12.3 152.0 
- MAy1 26:00.00  83:36.00 100 500 6.2 46.8 41.5 235.1 ~9.6 21.1 8 =-54.1 277.7 139.0 
= MA2-—s 2605.00 = 83:42.00 50 7S 1.1 23.9 -14.7 27.4 4.0 28.1 43.4 139.4 171.0 


~ Instrument and water depth in seters. 


- All speeds inaas . 


- Instrument depth on model morings = ‘s(200 m). 
- Model data has been truncated to match observed data. 


Ff 


INST. WATER 

D LAT LONG DEPTH «3=s«CDPEPTHCCU+EANCWH 
MW2 =_-25:24.00 + 85:24.00  =(172302Csi32K6—i(‘<ié‘é« = D=wCW 0.3 
Mjch «= «2536.03 =: 8529.97 1565 3200 1.0 39.9 
MA2 = =s-«:25:36.00 84:48.00 1172 2143 0.0 3.2 
MAM =—s«-25:42.90 84:53.10 1100 1697 -0.6 8.5 


- Instrument and water depth in meters. 


~ All speeds in one, 


gt 


- Instrument depth on sodel soorings = 200 + (lower layer depth). 
- Model data has been truncated to match observed data. 


simulated current meter data at approximately the same 
locations from the 0.1 degree port-forced experiment 201/13.0. 
Figure 40 shows the locations of the moorings; moving eastward 
from deep to shallow water the simulated data sets V, W, X, Y, 
Z correspond to moorings G, A, C, D (DA), F, respectively. 
Since the ocean model is a simulation, rather than a 
prediction, it is not possible to make a day by day comparison 
of the data sets. However, a realistic simulation should 
behave similarly to the actual Gulf and, therefore, it is 
appropriate to compare various statistics from the observations 
and the simulation. The major limitations of the comparison 
are that the ocean model does not have wind forcing, and the 
observational data covers less than two eddy cycles with gaps 
due to instrument failures. Another possible limitation of 
comparisons at a small number of locations is that relatively 
minor errors in the position of a simulated circulation feature 
might lead to large differences from the observations; however, 
that does not appear to be the case in this instance. 


Table 1 summarizes the statistics at the three moorings 
(C, DA, F) and three simulated locations (X, Y, 2) on the 
shelf, and at mooring E which is north of mooring C. In all 
cases the currents are significantly stronger in the simulation 
than in the observations. Figures 41 and 42 are a comparison 
of observed and simulated kinetic energy spectra at moorings C 
and F, but note that the use of a linear rather than the more 
usual logarithmic scale for variance tends to magnify the 
differences between the data sets. The simulation at C, MOXl, 
has a very strong signal at a period of about 400 days and very 
little energy at periods less than 30 days. A port-driven 
experiment should not have much energy at short time scales, 
Since wind forcing is dominant at these scales. But the energy 
variance peak indicates that the Loop Current, with an eddy 
shedding cycle of about 400 days, acts directly at MOX1. The 
observations at C, MOC2, have a comparatively level spectra 
with a minor peak at about 85 days and with energy at periods 
down to 5 days. The Loop Current clearly does not reach 
location C in the Gulf, the minor peak in the spectra might, 
however, be due to the passage of filaments that are often 
observed on the wall of the Loop Current. The simulation at F, 
MOZ1, is far too energetic and still shows a strong influence 
from the Loop Current. In fact, M021 is a better fit than MOX1l 
to MOC2 on the shelf edge. Data from the corresponding 
observational location, MOF1, produce a very flat spectrum and 
very little energy variance, as is expected in water that is 
only 50 m deep Overall the ocean model exhibits a very strong 
Loop Current effect over the shelf area, an effect that is not 
present in the observations. 


Table 2 summarizes the statistics at the two moorings (G, 


A) in deeper water and at the cooresponding simulated 
locations (V, W). Only deep currents are included, this is 
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unfortunate since the ocean model would probably show good 
agreement with the observations at about 100 m depth. The 
Simulated deep currents are significantly weaker than those 
observed in the Gulf, particularly at mooring G which is in 
3,000 m of water directly in the path of the Loop Current. The 
deep currents for experiment 201/17.0 are quite similar to 
those for 201/13.0, and Figures 29 to 34 indicate that the 
simulated currents are probably stronger between V and W than 
they are at these locations, and that the strongest simulated 
deep currents under the Loop Current are further south. 
However, even if it is assumed that the ocean model is 
producing the same deep water circulation features as those 
observed, but with a shift to the south, there would still be a 
factor of 2 to 3 difference in current speed ranges. The 
simulated deep currents are very sensitive to the amount of 
deep transport through the Yucatan Strait, but there is little 
data available on the vertical distribution of flow through the 
Strait. Many of the initial simulations of the Gulf had only 4 
Sv inflow in the deep water, more recent simulations use 10 

Sv. The data at moorings G and A suggest that an even larger 
deep transport might be appropriate. One of the main reasons 
for locating moorings at G and A, in deep water, was the need 
for data that could be used to tune the ocean model's deep flow 
characteristics. This is a good example of cooperation between 
the physical oceanography and ocean modeling programs at MMS. 
Simulations in Year 3 of the modeling program may use a larger 
deep inflow transport, however surface currents are more 
important that deep flow so these simulations will only be 
delivered to MMS if they give rise to more realistic upper 
layer currents. 


In Year 3 of the modeling study, now in progress, an 
effort will be made to develop an efficient layered ocean model 
that will allow layer interfaces to effectively intersect the 
bottom topograhy. It is likely that more than one workable 
solution to this problem can be found if computer cost is not 
taken into account. In fact, the major difficulty will be in 
finding a practical solution that does not impose too high an 
overhead for its use. When intersection occurs in a 
coventional layer model the layer thickness becomes negative 
which is clearly unphysical, leading to unrealistic results 
and, if the situation persists, catastrophic failure of the run 
due tc undamped instabilities. The obvious solution of setting 
a minimum layer thickness at or about zero does not work 
because, (a) it leads to loss of mass, and (b) clamping the 
layer thickiess induces dispersive ripples in the interface at 
the intersection point (i.e., we have an unresolved boundary 
layer). One promising approach to the layer intersection 
problem is to insure positive layer thicknesses via ‘Flux 
Corrected Transp~~t', a technique that was originally developed 


for fluid pre’ es with shocks (Book et. al., 1981). In this 
method the cor ::iwity equation is solved for layer thickness 
using two ¢ scent sets of transports, one obtained via a low 
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order (nighly dispersive) scheme guaranteed to give monotonic 
results and the other via a standard (ripple prone) high order 
scheme. The low order scheme used alone would prevent layer 
intersection, but it very rapidly damps out circulation 
features and therefore would not produce realistic 

simulations. Instead, the final layer thickness at each point 
is a Linear combination of the two solutions, chosen to be as 
close as possible to the high order solution. Away from areas 
of layer intersection the high order scheme will dominate, but 
near intersections just sufficient contribution from the low 
order scheme will be used to ensure overall positivity. In 
other words, bottom topography is still confined to the lowest 
layer, but that layer can get very thin so there is effectively 
no contribution from the deep layer over the shelf and no limit 
on how shallow the bottom t®pography can be. This method has 
been used with some success for interfaces that intersect the 
surface (Bleck et. al., 1983), but is usually implemented in 
fully explicit models. Explicit free surface ocean models are 
not used for long term simulations because they require very 
Short time-steps. To be practical, therefore, flux-corrected 
transport must be modified to work with the present 
semi-implicit time stepping scheme, or perhaps a split-explicit 
version of the model must be developed. 


In order to demonstate the viability of using flux- 
corrected transport for layer intersection problems, a two 
dimensional explicit version of the ocean model for flux- 
corrected transport was developed by JAYCOR for NORDA in the 
summer of 1985. The goal of the project was to run a numerical 
experiment that is simple to describe but very difficult to 
perform with any kind of conventional ocean model. Two basins 
are separated by a sloping wall, i.e. a seamount, and at time 
zero the first basin is partially full of water with the second 
empty. As the experiment progresses more water is added to the 
first basin until it overflows the seamount and the second 
basin starts to fill with water, eventually both basins are 
full of water. Conventional layered ocean models cannot even 
Simulate the first stage of the experiment, which requires the 
air-sea layer interface to intersect tne sloping seamont 
topography. No conventional ocean model can simulate the 
second stage of the experiment, in which the second basin fills 
with water. In order to run the simulation in a reasonable 
amount of time, the experiment was performed in a parameter 
space quite different from a typical ocean basin. For example, 
the water was only 100 m deep with a very gently sloping 
seamount. Moreover, the simulation only indicates that the 
method can allow layer interfaces to intersect topography to be 
freceswer flux-corrected transport must be extended to allow 

onger timesteps than those use in explicit models. 


Figures 43 to 51 show layer depth every 6 hours for the 


full 2 days of the experiment described above. At time zero 
the air-sea interface is defined across the entire region, but 
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the layer is very thin in the second basin. Water is added 
Slowly at the center of the first basin and the water level 
rises. At 12 hours (Pigure 45) the level is about 15 m above 
the top of the seamount, and water has started to run down its 
opposite side into the second basin. At later times the level 
in the first basin continues to rise because the water is being 
added more quickly than is allowed to run into the second 
basin, but the pressure head directly above the seamount 
continues to be about 15 m. At 42 hours (Figure 50) the water 
level in both basins is above the seamounts, but there is still 
a notch in the air-sea interface above the seamount. In Figure 
51 the levels in the two basins are almost equal. The reason 
for the fairly large pressure head above the seamount and the 
Slowness of the flow into the second basin is that flux 
corrected transport requires that velocities in a layer tend to 
zero as the layer thickness goes to zero. This is achieved by 
using layer thickness dependent bottom friction that goes to 
infinity as the layer thickness tends to zero. In practice 
this means that thin layers become thicker only slowly; 
however, this happens quickly enough for layer interface to 
move up (and down) the topography in tue first layer at a 
reasonable rate. The notch seen in the interface is Figure 50 
is an artifact of the flux corrected transport method, it 
illustrates that the method may not produce realistic results 
near the point at which the interface intersects topography. 
This will almost certainly be true for any method that allows 
layers to intersect topography. In the Gulf of Mexico, 
therefore, velocities will be questionable at the shelf break 
(particularly in the lowest layer) but realistic over the shelf 
and over deep water. 


20 


VII. 


10. 


ll. 


12. 


REFERENCES 


Bleck, R.C., Rooth, D.B., Boudra 1983. “Wind-Driven 
Spinup in Eddy-Resolving Ocean Models Formulated in 
Isopycnic and Isobaric Coordinates". Rosenstiel School 
of Marine and Atmos. Sci., U of Miami, Florida. 


Book, D.L., J.P. Boris, S.T. Zalesak 1981: 
“Plux-Corrected Transport® in D.L. Book (ed), Finite 
Difference Techniques for Vectorized Fluid Dynamics 
Calculations, Springer-Verlag. 


Brooks, I.H. 1979. Fluctuations in the Transport of the 
Florida Current at Periods Between Tidal and Two Weeks. 
J. Phys. Oceanogr. 9:1048-1053. 


Elliott, B.A. 1979. Anticyclonic Rings and the 
Energetics of the Circulation in the Gulf of Mexico. 
Ph.D. Dissertation, Texas A&M University, 188 p. 


Hellerman, A. 1967. An Updated Estimate of the Wind 
Stress of the World Ocean. Mon. Wed. Rev. 95:607-628. 


Hurlburt, H.E. and Thompson, J.D. 1980. A Numerical 
Study of Loop Current Intrusions and Eddy Shedding. J. 
Phys. Oceanogr. 10:1611-1651. 


Leipper, D.F. 1970. A Sequence of Current Patterns in 
the Gulf of Mexico. J. Geophys. Res. 75:637-657. 


Marsh, J.G., Cheney, R.E., McCarthy, J.J. and T.V. 
Martin. 1984. Regional Mean Sea Surfaces Based on 
oo and SEASAT Altimeter Data. Marine Geodesy 

4 5-402. 


Maul, G.A. and Herman, A. 1984. Mean Dynamic Topography 
of the Gulf of Mexico with Application to Satellite 
Altimetry. Marine Geodesy (to appear). 


Merrell, W.J. and Morrison, J.M. 1981. On the 
Circulation of the Western Gulf of Mexico, with 
Observations from April 1978. J. Geophys. Res. 
86:4181-4185. 


Niiler, P.P. and Richardson, W.S. 1973. Seasonal 
Variability of the Florida Current. J. Mar. Res. 
31:144-167. 


Rhodes, R.C., Thompson, J.D. and Wallcraft, A.J. 1986. 
The Navy Corrected Geostrophic Wind Data Set for the Gulf 
of Mexico. NORDA tech. rep. (to appear). 


21 


REFERENCES (Cont.) 


13. 


14, 


15. 


16. 


17. 


18. 


19, 


Schott, F.A. and Zantopp, R.J. 1985. On the Seasonal and 
Interannual Variability of the Florida Current. Science 
227:308-311. 


Sturges, W. and Blaha, J.P. 1976. A Western Boundary 
Current in the Gulf of Mexico. Science 192:367-369. 


Thompson, J.D. 1986. Altimeter Data and Geoid Error in 
Mesoscale Ocean Prediction: Some Results from a Primitive 
Equation Model. J. Geophys. Res. 91:2401-2417. 


Wyckhouse, R.J. 1973. Synthetic Bathymetric Profiling 
System (SYNBAPS). Naval Oceangraphic Office Tech. Rep. 
TR-233. 


Vukovitch, F.M. and Maul, G.A. 1984. Cyclonic Eddies in 
the Eastern Gulf of Mexico. J. Phys. Oceanogr. (in press). 


Waddell, E. 1986. Gulf of Mexico Physical Oceanography 
Program Final Report: Years 1 and 2. (Final Report in two 
volumes by SAIC submitted to Minerals Management Service, 
Metairie, LA. Contract No. 14-12-0001-29158). 


Wallcraft, A.J. 1985. Gulf of Mexico Circulation 
Modeling Study Annual Progress Report: Year 1. (Progress 
report by JAYCOR submitted to Minerals Management Service, 
Metairie, LA. Contract No. 14-12-0001-30073). 


22 


FIGURES 


23 


FIGURE 1: (a) Instantaneous view of the interface deviation 
in a two-layer simulation of the Gulf of Mexico driven from 
rest to statistical equilibrium solely by inflow through the 
Yucatan Straits (Experiment 9). The contour interval is 25 
m, with solid contours representing downward deviations. 

(b) Depth of the 22 degree isothermal surface, 4-18 August 
1966 (Alaminos cruise 66-A-i1), from Leipper (1970). The 
contour interval is 25-m. 


24 


FIGURE 1 


(a) 
¥" “. . \ as 
1 \ Ss 
~ i ‘= 
; i ‘" “ 
| \ — sf 
(cr) | JE 
~.! 
\ 
SZ 
\ —_ 
== \ : \ Se 
= \ 7 
0 ‘ ee _ : 7 
0 (an) 1750 
(b) 


25 


BEST COPY AVAILABLE 


FIGURE 2: (a) Interface deviation from the Gulf of 
Mexico simulation at model day 1970 after an eddy has 
separated from the Loop Current and propagated westward. 
(b) Ninety days later the major anticyclonic eddy has 
developed into a counter-rotating vortex pair in the 
western Gulf. The cyclonic vortex is to the north and the 
anticyclonic to the south. 
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FIGURE 3: Counter-rotating vortex pair in the western Gulf 
of Mexico as shown by the depth of the 15 degree isotherm 

(in meters), observed in April 1978. The cyclonic vortex is 
to the north and the anticyclonic to the south (from Merrell 


and Morrison, 1981). 
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FIGURE 4: Sea surface height variability for the Gulf of 
Mexico. (a) Based on about 16,000 GEOS-3 and SEASAT 
cross-overs, spanning nearly four years (from Marsh et. 
al., 1984). (b) Based on an ocean model simulation with 
port forcing only (Experiment 40), measured over three 
eddy cycles at statistical equilibrium with the free 
surface sampled every ten days for a total of over 300,000 
"“observations”. 


FIGURE 5: Sea surface height (a) variability and (b) 
mean, for the Gulf of Mexico. Based on all available 
hydrographic, STD and XBT data at over 16,000 stations, 
with substantial filtering (from Maul and Herman, 1984). 


FIGURE 6: Mean sea surface height for the Gulf of 

Mexico. Based on an ocean model simulation with port 
forcing only (Experiment 40), measured over three eddy 

Ne at statistical equilibrium. The contour interval 
s cm. 
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FIGURE 7: Instantaneous view of upper layer averaged 
velocities from Experiment 68 on model day 3918. Vectors 
are only plotted at every second model grid point, 1.@., 
every O. 4 degrees. 


FIGURE s : Instantaneous view of upper layer averaged 
velocities from Experiment 68 on model day 4038. Vectors 
are only plotted at every second model grid point, i.e., 
every 0.4 degrees. San 


FIGURE 9: Instantaneous view of upper layer averaged 
velocities from Experiment 68 on model day 4158. Vectors 
are only plotted at every second model grid point, i.e., 
every 0.4 degrees. ae 


FIGURE 10: Instantaneous view of upper layer averaged 
velocities from Experiment 68 on model day 4218. Vectors 
are only plotted at every second model grid point, i.e., 
every 0.4 degrees. ie 
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FIGURE 8 
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FIGURE 9 
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FIGURE 10 
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FIGURE 11: Seasonal climatology for winter (Dec., Jan., 
Feb.) from ship observations (Elliott, 1979) and’from the 
Navy Corrected Geostrophic Wind data set (1967-1982): (a) 
wind stress 1967-1982, (b) wind stress curl 1967-1982, (c) 
ship wind stress, and (d) ship wind stress curl. ~ 


FIGURE 12: Seasonal climatology for spring (March, April, 
May) from ship observations (Elliott, 1979) and from the 
Navy Corrected Geostrophic Wind data set (1967-1982): (a) 
wind stress 1967-1982, (b) wind stress curl 1967-1982, (c) 
ship wind stress, and (d) ship wind stress curi. ~ 


FIGURE 13: Seasonal climatology for summer (June, July, 
Aug.) from ship observations (Elliott, 1979) and from the 
Navy Corrected Geostrophic Wind data set (1967-1982): (a) 
wind stress 1967-1982, (b) wind stress curl 1967-1982, (c) 
ship wind stress,'and (d) ship wind stress curl. - 


FIGURE 14: Seasonal climatology for fall (Sept., Oct., 
Nov.) from ship observations (Elliott, 1979) and from the 
Navy Corrected Geostrophic Wind data set (1967-1982): (a) 
wind stress 1967-1982, (b) wind stress curl 1967-1982, (c) 
ship wind stress, and (d) ship wind stress curl. 
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FIGURE 15: Instantaneous wind stress and wind stress curl 
from the Navy Corrected Geostrophic Wind data set, for 0000 
and 1200 GMT on 14 January and 0000 GMT on 15 January 1976. 


FIGURE 16: Instantaneous wind stress and wind stress curl 


from the-Navy Corrected Geostrophic Wind data set, for 0000 
and 1200 GMT on 14 July and 0000 GMT on 15 July 1976. 
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FIGURE 17: Instantaneous view of the interface deviation in 
the winter from two simulations driven solely by clima- 
tological winds. (a) Experiment 31 driven by seasonal ship 
winds, and (b) Experiment 202/13.0° driven by monthly winds 
from the Navy Corrected Geostrophic Winds data set. The 
contour interval is 12.5m, with solid contours representing 
downward deviations. 


FIGURE 18: Instantaneous view of the interface deviation in 
the spring from two simulations driven solely by clima- 
tological winds. (a) Experiment 31 driven by seasonal ship 
winds, and (b) Experiment 202/13.0°driven by monthly winds 
from the Navy Corrected Geostrophic Winds data set. The 
contour interval is 12.5m, with solid contours representing 
downward deviations. ~ 


FIGURE 19: Instantaneous view of the interface deviation in 
the summer from two simulations driven solely by clima- 
tological winds. (a) Experiment 31 driven by seasonal ship 
winds, and (b) Experiment 202/13.0°driven by monthly winds 
from the Navy Corrected Geostrophic Winds data set. The 
contour interval is 12.5m, with solid contours representing 
downward deviations. 


FIGURE 20: Instantaneous view of the interface deviation in 
the fall from two simulations driven solely by clima- 
tological winds. (a) Experiment 31 driven by seasonal ship 
winds, and (b) Experiment 202/13.0° driven by monthly winds 
from the Navy Corrected Geostrophic Winds data set. The 
contour interval is 12.5m, with solid contours representing 
downward deviations. 
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FIGURE 21: Instantaneous view of the interface deviation in 
the winter from Experiment 202/13.0 driven solely by winds 
averaged by month from the Navy Corrected Geostrophic Winds 
data set. (a) winter of 1970/1971, and (b) winter of 
1971/1972. The contour interval is 12.5m, with solid 
contours representing downward deviations. 


FIGURE 22: Instantaneous view of the interface deviation in 
the spring from Experiment 202/13.0 driven solely by winds 
averaged by month from the Navy Corrected Geostrophic Winds 
data set. (a) spring of 1971, and (b) spring of 1972. The 
contours interval is 12.5m, with solid contours representing 
downward deviations. 


FIGURE 23: Instantaneous view of the interface deviation in 
the summer from Experiment 202/13.0 driven solely by winds 
averaged by month from the Navy Corrected Geostrophic Winds 
data set. (a) summer of 1971, and (b) summer of 1972. The 
contour interval is 12.5m, with solid contours representing 
downward deviations. ~ 


FIGURE 24: Instantaneous view of the interface deviation in 
the fall from Experiment 202/13.0 driven solely by winds 
averaged by month from the Navy Corrected Geostrophic Winds 
data set. (a) fall of 1971, and (b) fall of 1972. The 
contour interval is 12.5m, with solid contours representing 
downward deviations. ~™ 


D2 


PIGURE 21 


INTERFACE DEVIATION 
DATE = 3586/1970 DH = 125(M) LAYER =1 G. OF MEXICO 20232:2: 11.0 


SIN ' i ' ' La ' ' ' | Lj ' 1 ' q T Zs 1 


188 i —— il l j ee — 
oer MIN = ~—4.33E+01 MAX = 1.11E+02 sow 
NwOPD4 323 02/13/86 
INTERFACE DEVIs TION 
DATE = 353/1971 DH = 12.5 (M) LAYER = 1 G. OF MEXICO 20232:2: 11.0 
JIN - a ae ae T Tt a a eee ee ee a | 
16N 1 i — a baal i i 
96W MIN = -§.10E+01 MAX = 1.13E+02 sow 


53 NOPD4 323 02/13/86 


GEST COPY AVA ARTE 


oil WoL 


31N . 


INTERFACE DEVIATION 
DATE = 083/1971 DH = 125 (M) 


FIGURE 22 


LAYER = 1 G. OF MEXICO 20232:2: 11.0 


— a ’ ss qT 


' Li “. T ' ' i i 7 ' 


+ 


p 
i l i 1 j ij i 1 l i 
ia) MIN = ~4.84E+01 MAX = 9.42E+0i bow 
wOPOA 323 02713786 
INTERFACE DEVIATION 
DATE = 078/1972 DH = 12.5 (M) LAYER = 1 G. OF MEXICO 20232:2: 11.0 
1 1 | ' 1 | i 1, 17 1 Li 1 Li v i 7, 1 


1 


L i j lj in 


i j L 
MIN = ~5.37E+01 


j j 
MAX = 1.05E+02 aow 


NOPOA 323 02713786 


BEST COPY AVAILABLE 


D4 


INTERFACE DEVIATION 
DATE = 173/1971 DH = 12.5 (M) 


FIGURE 23 


LAYER = 1 


G. OF MEXICO 20232:2: 11.0 


31N a ao T T 


a 


i ' ' 1 


j 1 


18N i 


=e 


i. 


= 


' ' Tt ' 1 


i 1 L 1 l 


INTERFACE DEVIATION 


DATE = 168/1972 


DH = 12.5 (M) 


MIN = =3.13E+0) 


i 
MAX = 1.15E 


LAYER = 1 


+02 


NOPOA 323 02/13/86 


G. OF MEXICO 20232:2: 11.0 


1 


T 


1 


, 


T 


1 


q 


L 1 l a. 1 m | a. 1 


5.02E+01 


16N — 


an 
96¥ = 


j 7 
MAX = 9.41E+01 


y 


MIN Bow 


BEST COPY AVAILABLE 


NORPDA 323 02/13/86 


FIGURE 24 
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FIGURE 25: Bottom topography and coastline geometry for 
Gulf of Mexico model on 0.1 degree grid. The contour 
interval is 500 m, and there is a flat shelf between the 500 
m isobath and the coast. The section of the Caribbean shown 
is treated as land by the model, the position of the inflow 
port is marked by the termination of contour lines in the 
Yucatan Strait. 
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FIGURE 26: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/13.0 on model days (a) 3150, 
and (b) 3240. Since this experiment has no wind forcing the 
assigned dates (249/1968 and 339/1968) are arbitrary. 
Vectors are only plottéd at every’ sécond model grid point, 
i.e., every 0.2 degrees, and all velocities greater than 50 
cm/sec are plotted as 50 cm/sec. 


FIGURE 27: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/13.0 on model days (a) 3330, 
and (b) 3420. Since this experiment has no wind forcing the 
assigned dates (063/1969 and 153/1969) are arbitrary. 
Vectors are only plotted at every second model grid point, 
i.e., every 0.2 degrees, and all velocities greater than 50 
cm/sec are plotted as 50 cm/sec. 


FIGURE 28: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/13.0 on model days (a) 3510, 
and (b) 3600. Since this experiment has no wind forcing’ the 
assigned dates (243/1969 and 333/1969) are arbitrary. 
Vectors are only plotted at every sécond model grid point, 
i.e., every 0.2 degrees, and all velocities greater than 50 
cm/sec are plotted as 50 cm/sec. 
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FIGURE 29: Instantaneous view of the (a) upper and (b) 
lower layer averaged velocities from Experiment 201/17.0 on 
model day 2070. Since this experiment has no wind’ forcing 
the assigned date (264/1965) is arbitrary. Vectors are only 
plotted at every second’ model grid point, i.e., every 0.2 
degrees, and all velocities greater than 50 cm/sec are ~ 
plotted as 50 cm/sec. 


FIGURE 30: Instantaneous view of the (a) upper and (b) 
lower layer averaged velocities from Experiment 201/17.0 on 
model day 2160. Since this experiment has no wind forcing 
the assigned date (354/1965) is arbitrary. Vectors are only 
plotted at every second’ model grid point, i.e., every 0.2 
degrees, and all velocities greater than 50 cm/sec are 
plotted as 50 cm/sec. 


FIGURE 31: Instantaneous view of the (a) upper and (b) 
lower layer averaged velocities from Experiment 201/17.0 on 
model day 2250. Since this experiment has no wind’ forcing 
the assigned date (079/1966) is arbitrary. Vectors are only 
plotted at every second’ model grid point,’ i.e., every 0.2 
degrees, and all velocities greater than 50° cm/sec are 
plotted as 50 cm/sec. 


FIGURE 32: Instantaneous view of the (a) upper and (b) 
lower layer averaged velocities from Experiment 201/17.0 on 
model day 2340. Since this experiment has no wind’ forcing 
the assigned date (169/1966) is arbitrary. Vectors are only 
plotted at every secOnd model grid point, i.e., every 0.2 
degrees, and all velocities greater than 50 °cm/sec are ~ 
plotted as 50 cm/sec. 


FIGURE 33: Instantaneous view of the (a) upper and (b) 
lower layer averaged velocities from Experiment 201/17.0 on 
model day 2430. Since this experiment has no wind forcing 
the assigned date (259/1966) is arbitrary. Vectors are only 
plotted at every second’ model grid point, ’'i.e., every 0.2 
degrees, and all velocities greater than 50 cm/sec are 
plotted as 50 cm/sec. 


FIGURE 34: Instantaneous view of the (a) upper and (b) 
lower layer averaged velocities from Experiment 201/17.0 on 
model day 2520. Since this experiment has no wind forcing 
the assigned date (349/1966) is arbitrary. Vectors are only 
plotted at every second’ model grid point, i.e., every 0.2 
degrees, and all velocities greater than 50 cm/sec are 
plotted as 50 cm/sec. 


63 


FIGURE 29 


DATE = 264/1965 
G. OF MEXICO 20132:2: 17.0 


| oS u/s 

jain T T T a oe pony T T T T T T 

{ - + + 

{ 

{ 

t 

i - + . 

' =. 
bie ener a4 

i we nase 

‘ ey a A, 

oa ate 

‘ 

‘ 

i 

‘ 


Fa 8 
™= 


aad 


4420) 
7 AAAP A> DOP 
Sd 7 4 


4s 
Dd 


aN 


1 l 


‘LOTTED SPEED = 1.04 (M/SEC) 


NOPDA 323 02/28/86 


DATE = 264/1965 
G OF MEXICO 20132:2: 17.0 


LAYER 2 CURRENTS 


— 


7 
Sd v 
<<er 44na4 4 


s -=4+> 
5 797 7p 
“ oY eave A Mm 
Rowe V4 >>ar oe << <9 44440 © ORR 44 
as rere ert é 


> a> 
<tr TATA A 
> Paryars ATATe > aad EPR ay Pheer 4 apc 4 
eng aae Arann? > eaten Ae +¢ _ 
06 <> Pd Oe eee4 44> FPP? babe 44 Bae F 224 Vice 
a¥ ape gy > ive Ses? eaa?* paiVve <n AY 
avre Vos bead a a ON ee bbe & bee ay 
oo rev 4 


» esac 4 4 4 
oe 
< <erer? + 4aq 4% 


> at 
a tie 


7 

Ale oo 2 
te Oe eee < 
© <4 eerenn t 


a 

a” 
‘ 

v «ser 


Yq? 


& 


4 4 cs 

> ° < . 
yola4 pa 4940) Casta SvY bepneP 
? 2 74 %re<e Verto 
torte A Chae ae i A oe 

YP Anbbel ye ARV de 
$ ev > 
7. 


> <Vv yew ¥ &-m> >>> 
ons ech ava 
.——ooe 


a 


aaa 


4 Serttitt 


aa 
r ee Ap wadsosseny 
on pais <a 
he deh werent. 2 
4? @>>* 
+eeeeer aq*¥ a 
eer 446444 + ha rer VV ean ReneT< 
qeVece VAN be Keeg ee Sree ve > tr TAH4 D> TP? & 
TTP © whi pee ge KY PRE HKe ryrvaat mae 
VSPEPPRR Fe dy oe <i « <<< A 
qq 44444 IUk aaeee © F< 44 © open chPan 
«co 


Ex rays 74 


Gd, aaaee? 7° € 44 4 4 CECTTT RATA 


<4444¢¢ 
i 


+ 
7 


ae 
4 4ee4¢ 


wv" 


+ 
dl 


77777 
< 


> 
< 


TASS 
a ete 


Ree 
.2it 


o 


Mii 


a® 


A 


j 


— 


/ j 


— 
MAX PLOTTED SPEED 


+ fa 


j 1 
= 0.23 (M/SEC) Bow 


NOPOA 323 02/28/86 


64, 


LAYER 1 CURRENTS 


FIGURE 30 


DATE = 354/1965 


os u/s G. OF MEXICO 20132:2: 17.0 
3iN T Y T T T T T =F a T T T T T T 


we 
“4 
(7; 


aa 


ZR 
‘és ay 


* 


Views acarrr> 
re ‘ 
«a? 


stteecace? 
~ oseee< 
vi¥ardo2v 


: i eatert : 
piri 
eg Yes l 2... d 
? 
16N | j j | ee | i L 1 j] i 
o6W MAX PLOTTED SPEED = 1.63 (M/SEC) eow 


2 CURRENTS 


NOPOA 323 02/28/86 


DATE = 354/1965 


05 u/s G. OF MEXICO 20132:2: 17.0 
3IN T T T T T T ’ aay a | a T T 7 
= > + > + + + + 7 7>>>? * <c<> + 7 * = 
Ay eye RB aN a & 
> 233 ve 
7 > > swe f. vetteares tg * SS" 
Sg . 
— . . ger aaasi S5DYSD >> fm Tes + eee cee + + om 
5 4 4 ae 
4494993488099 9 9 35 Ct S a Nasnnas SV en A Ap 
> a +4 
Ae Ae Le ee ee a pa4ne 9% aga> Vtaa> & ® @enae 4 P44 § 
= + ine > ZAREREeVV YE 7 Aa 4A > a 74 + = 
ffege>? 2O9 << +> Pee? Pr PF F FHA a>? > > >> 4 4 gad SA A 
pits: "aay pa Vaece 7 ee Pz Eeeme MP> coasaer> a 
iy a >pa4> 72> berrr > we Se, Pree rae e444 6% Conan 474% 4 +99 ‘ r » * 
errs siashiereetneaiteg Lat’ Sia Bee al seta ab fas hee tps hs 
-— Paar aa rpaaaee® AA 77744 A ‘ aes 4 ° ° 
fies te a tase ve «293445 S* se TRIS 5593> 4 444 
ate + 00 ae a » ° <<<? , >.> ates werrrrven” oa Gee . z, A RA 
¥ > es Bars see Eos rar > 7 eee . 
win > a a >eneeeee «3 wyis hy ae 
49pavseee ee en , avvNSa eee sere annasehe” vie + 
ME AT Td Ahbbbdooeeeccae € F € soe i" s¥93, a4 vt A 
«aa teen Aan 
< eae ae y aaeaeet esas ake w*gdss vv Sent 
° < 
“6+ ~@as32> -“* egcervos, aie, oT lel 
ae ehedddda *¥ » M4 anee 
errhp Veecee oat "Fan se <4 roa 
weweenhaes he hae ‘9 Le | 
ae my _piebes eee Nytsze end >» 
‘ aes sot th ah al wR 704s N 
‘y + e+ AAT CA A a 
<~ z n be erp e Revd ® arcvae 
aé«é ? _ ia i " eases Ve atp anes cteer> Vaeeetan 
* > nce bey - Sorry 
er 4 pew nn #, eeeennera ne a "tisrenpés 
> wees 7A ep Pe 
rs 466 < ey ee < opty" a "e >> vane pit bh. Zs 
<“te?¢ * - 
3 vy VYsse AF bs 67 ra r, Nad A AAADP?? * res f 
>’ see perernas7. a* “945 
Py a> ata aieie” <* piss > aver "haan 
< 41> +, 4+ VA - > PSO >> «e ¢ aa 
er wet &< <“¥y vec <« <@ 27 
+ . vy ve a 
*“ry <“ts RRR A > > ~“ 
me vw ee" 4a" A 
ard ** sevvnedie a? 4 } 
* ein, ee eheens eee aaz 
. Tes 4 cme J s+ “vy + + + -4 
Trev beaver eeteeew' VY", 
Rye we £6 be nvr ng 
fire yet hv iwsowe 4 
4 wybeceten 
- + eVFCe aaVVe «4 bee t & 7 + + 
et 4 rF xy 
aac 444s ? 
<<< 
1@N 1 rT j L iL lj 1 i i | 1 l lj jl = 5 
968 MAX PLOTTED SPEED = 0.27 (M/SEC) sow 


BEST COPY AVAILABLE 


wOPDA 323 02/27°786 


65 


1 


FIGURE 31 


DATE = 079/ 1966 
G. OF MEXICO 20132:2: 17.0 


' ' qT Li 


7 


<<<? 


4ve 
3 


° 
> 


iy,’ 
ee 
Pe << 


? 


16N , i 1 l j i lj i i i l l 
96W MAX PLOTTED SPEED = 1.05 (M/SEC 60" 
NORPDA 323 02/28/06 
LAYER 2 CURRENTS DATE - 079/1966 
05 u/s G. OF MEXICO 201232:2: 17.0 
3iN T T T T T T T T T T pee pee 2 T T = 
— > + + + . + + + pias, £4 cen «+ + 
hte” > >» > ov 6 
77? ~~. a 7 > Pepeee ¥ << «<< 
-€ q ve Pmvy < e< e 
= + + 794 art 4 Vetoes o> > 7 wv aay SV Se n> ee te 
as i 44449 *+ +4 v yw ¥en ae ie 
“77>? 7 ae* ¥ 4 SF + oe aeeet tap teay S"4 a a a daa gant, § 
>4 72? a*a ah 004 vv 
r>PreT > PF Petar H FY F FF4AT FT FePae F*? 2B * 244 FF va 
— +J/4 Tr>23 tA Pe > < 
g ° U2 per hes cae vee PO DPT DST Sane 2 eee 24 WY & Sup 
7 ' * 
A 4 44444 Se o> &® DOTTET? FOF a tae de ee 
phen hea 2 > > g > 7 > 744 ae te 
AB ehenegq> owes F>>09 77> > + + 7PSr> [At777 2 > PP Pad sae wean 
pram +4 Ree PTD v7 days prees*sa >» 4 ase 6 
RR AAD Y Ge Rabe reer OMe o4e FF > Fo att < APrad ttegeae a A 
7 “aclh, 7 > Socdbant My ?. > 8S « >Paaare es 4 > see 4 
A - 77 > ® 
- “14 PF S or¢ Vasa qe<cevved <oeess vo" ’ 
an an Pos teepecitas a mint a «OA 4 ORE LK AED oe ere 
fet ectaaavel tee Ve 
A nnd peeceeenees << evVVvVdNe” pemee € € 4 4 44 4 he a’ wethe, . 
“dass cee cdas abe <« eebadevbe ee ats A AA 
4 > a ac « A 
# age eee Lin arse uae vk* - — ss 6a A, tote N . 
4 > tate oo peeeet? VY © eee nst F > eA 4he 4g FF 5 om > 
2 av ghee + . é feng a’ qnrm A eb a? FAS <P 
4 544 Anv tes See cee Pate e > 2 4 Venn nage 
rast BEE ARE tea tty ud Wes these atm eeeommcene 
44)4> «4 vf + *. neg 4a wav 44 vee bbs 4 14 8 ots _s 
4% a‘ ‘ <n ov we > > be pe Coes + @ >. oeee 
+3 72 vs v cr & ot VWteae 4dr of *s AMADA A COC Ee 
i Shh 2 4 24> <ivs 22° 
#48 phe woh RA ARR whe ag Whe 
sata 4 ¥ * & aavv « 
a4 4nav vt an ot ae Fe le 
peitay i afc) Vue 
ge’ by vv'l e ree tag Df haa 
SPorw 4 vf < + eeeee wiv 174 
tea, +. wore wenn 
Oe ina oH 
Weeeore }44 
Nesters ef! ¢ « 
Sheyet} 
‘ 
atvy 2 ? 
t a ti-7 
a + malar’ & + ~ 
wetae = 
.. 
\ * +4 
. “Fa, 
a ry + 7 
16N i I i 7% - j in i i i | = i i i 
o6W MAX PLOTTED SPEED = 0.21 (M/SEC) aow 


BEST COPY AVAILABLE 


wOPOA 323 02/28/86 


66 


FIGURE 32 
LAYER 1 CURRENTS DATE = 169/1966 


os u/s G. OF MEXICO 20132:2: 17.0 
3IN t 1 |  H 7 Tt 1 1 | i , T T tT t 1 , 
— + . ~ 
_ + 


whaad 
OTT prermes s, 


hw Nee, 


NN: 


4 ‘2 te, 
i el dl A 7, 
einninnt at DA hy 


? 


ae ee l l i 
Ww MAX PLOTTED SPEED = 1.15 (M/SEC) eow 


mOROA 323 02/28/86 


LAYER 2 CUS.RENTS DATE = 169/1966 


05 u/s G. OF MEXICO 20132:2: 17.0 
—_——_ oO 
3IN 1 | 1 | if ? 1 ia 1 i 1 T 1 | 7 i qT 
ma + a > a > > > ¢ 7? ~~ . ~ 
ayy??? 44 8 
a4 a 44 spe 
IIPIPD | y's ie AT? ® pas’ a7 Pe Bem 
A RRAA RD hd Vv mA4q 
~ a o er ' sehen IF > hg 23t ¢ Mp o> > TT = 
64s “yy > oa 0 ets . 
Fm a Tapectertan > 475 9 b> asaaaet SYS 4 2 pages coo 4 
Fangs ‘a> > 4n¢ <q7 > «<* 
7 > &@$4 4 beep FF FF FFF hel ea 7 ea — > we. vea aad v 
Te Veo eters BP Pate ase ee > 2 p52 Ps 4 Yana @_¥ 
4 -«*VVVY Arar VV@ var Savesves "sadvveed av 
4 4% @> bg = gs 49 32249 > ” "Ss Yoo AT 22 seen Tee eerne 4 ay th “ 
Ny Se Rees > . 
a o <= ane rere > 2 ae a re rey Sd we pon Sha a wee Bae a taney ¥ 4 
< 74 A a* Vara aag EPeyhh | oe ae ny 
heve Y &o>, oY faa, Iba,, 4 ww x2: 4 « "EY ew saat 4 aN 
Arann 4 és cow + ecervvevsetve © 4a? 
A perce a & * $ > 2a 2h, Se £44) 6 ee pM fad Wwe, « “47 7% 4 4 
_ nec lO sresqress spi A> av <ivrvev v T, Vanee 
A atVananae ee < <<@a54 ** + ate v weer? ¢ af & &W 4¥en wier’s - <** © 
a7 r <tredn 9 Vb da Aan "ESS, quvé v 
*aReT ibys <6 SS ces ella oy SG tee ew pyere < <6 
5 A p22 4Arege we agi Rs 7 < 
: d Stangee © << < Cerupeeeees < < CODD ¢ 4 ce Sats + ov 
ates” woececeoe Pte, eae Sa eaand 
“44 %ee ave vw & . ace ¢€ 4 4 4 48° © fH ae 7 a5** => a” » « ga a’ FP vanr>>> | 
- kn, . > is ~'s A> oee 2 he 
Bats sv pee Sth reese « + a4¢4 - me 5s wwe" 4 >. 
>a Fewey ty Ss Vienne 4 9 HS PL at mths 

aA ow “*reeree rie saan ha asanes ager < 4 94 ca by N et 
a7 * ew vee a AIF pin Bnee 28 Hen hs ys Pr POST n> > . 

° Trveeesee ’ ‘o9 7744? Sine” aa fA ey ¥ eas 4 Var an Ne ET. AAAA 9 FORD HEH 
tur, At lewer ~ enemy ¥ A 129 7 AP aay >>> 
tata ne Siaeeets 4499 en 1 ne 4 aaanzeea 4 ww 274 4 te 

. Ay 2eSeee >a + Aka “ Ane <* om tq *.4 
nae oe Saget 274 T5557 E Ms t Aes Ane © ee AY Vet vials 
ee aut oe & CERqR AAS ae £4 4V 

Ph Ban OO Venue » venenne wasttereh he ceeee <<) > « 

a® z - wee -- we abhue, 4" 

ge% Qr"t "a Soeew eae ee pe <x e << 0 seses 

ence «¥ -e7r OA «cacceee aan 
Re} hee ge OOF > 1" +s * netveg 
Rak Detes a - ~ 
247 SI a A + > 
as vec + > “es tt. vy 1 
; ¥ <b oa <y***3} > 
7 Sesacae -« «*® ers 4 ‘4 
a >>> & and die *, y 7asys 
= > - <nee> nee ~« a* + + * 
grrrw beavvave a ¥ 
«* ¢ £6$@& 44444 are « v4 
j++ <> staat en ‘ o 
Newt « Ly J 
mo +> ave hope bd Ny «ses, > > . 
«4 ée« bat be) 
thet bas ? 
«< 
18N 3 i ji i l i i j i i ra j i i 
oe MAX PLOTTED SPEED = 0.21 (M/SEC) 


wOPO* 323 02/28/86 


67 


BEST COPY AVAILABLE 


FIGURE 33 


LAYER 


oS u/s 


l 


3iN 


DATE = 259/1966 
G. OF MEXICO 20132:2: 17.0 


t ; 
\ Fie 
tye: 


44, 
+, Tm, 


~~ 


8 


“ha 


16N 


? 


= 


' ' 1 7 | 


L 


CURRENTS 


MAX PLOTTED SPEED = 1.18 (M/SEC) 


wOPO4 323 02/28/86 


DATE = 259/1966 


05 u/s G. OF MEXICO 20132:2: 17.0 
JIN T T T T T T t T T T aay T T T 
a * > + + . + > > w?>>? ote + > . ~ 
+ rr © & a, 
rs Tveereee > 
7? >» = 7? Peeaqg 4 > > Yee 
a 7. - <-* © 
+ + 440aga4a0>d ~ TT Sa TM hee ety 2 > 2 4 94 eT e . + _ 
rnd aye ee pag Ad a ad sel RA 
7, . ~~ a>7? &® Be 215 pervvaaat? aaa bin yy 4 apettsy, 
a, 4° €4V4, Ate “veces Vee 
23 4 Ries >? > aa” » “een WP a >» we Ye ES a 
— + rey << 17. eat PP4ashna vv « + + + 
ete , et Pe ies eer a Aa,“ cf a > “st v been 
<< <+ <eas ee ah aA 
a5 ry a~ be eceg sees Testis enn es resties 424 
a : A ad a LaPeer y pee tee ety yonrs ty "have t, i332 hae 
PY a @ ns a? + + 
= shiver: ioe yy ee aS a> seis 2 vs ra tan 
» AAds >» aver * Fae vr a vy mr Svaeate “an 
An > 446 —> 4a) KEM Raa > vy ae tt vee 
A 2» fr «4 Lorena dh? caVVVY 4 céaseed << Ph ecivvrse & & & a + <ve - orn 
= ap &® 3" AR eV Peaev bn Feteceeee piv + 
a * We 4 > FETT HP? PPA» P Peay + ee Tne $ rs ai ae S Pebee 
Aves NZ, ba Aa??? $5 * Kha 
n<ee anv FPPeEPAA AAAA TS reeves oe" “ qt baadananneeed cnt | <pre 
ares at 5 eh ee tong <« n°, 5 A 
 haae? .<3% o << ¢ ¢ eneee © cates secce vr¢ we ee - - 
oot seeoees © Ore bY Nes 2 oe - 
wehgis cee ee 444 Sesesstye ee) bones 242% 4“ 54 SR Saas | 
° 444ee suerte , Dm VAGead nh ae 


4 aeeeeeer err «© «4 4 he or tw 


Tt neh, aoe eerere “Ww? ee z 
anc Yeeue V¥Y & > SPRhnae © > Sikres ww 
V PAAR BR Or ey vbeeent 44 4 Iette aor? & A 
$34. Ge © any vt ww vitbe esses eee See 433" kel 
aim, -e 
eras, ee by VY 4 BO wetev * heat av 4 
44% . as te “s ‘ a 
atte vv vy Ve Fearne Seton nwt « }?? «4 => 
wher > eo ey Tistyr ary > 
mata YY 44 9 9 Pe 4 eC 66 bObmaan oS 
amu 4 
Tht S 44> 0 0's shoyu le eae 4 ¢i< «we 
> > aac 
SUS Yes 4b vo sheet By v'e ge 
—_ fray Vo aos wwitts Sew vv - 
: r A Joa 2% vy J Od v 7 . 
fv_a , A» eWreery WRARe © 
vw ave 4 <7 wh we G 
i > 4 + are «4 ch” garnet & 
> a «< ane 7 + > « 
pua* ¥ ¥ cass aos <#s 
RAR 66 KV ORG a 
Tr re yet bekiwee & & & 
ec 4 eh ve 
— + Weeeeacabe 4 sta, + + + - 
© <446%s ? 
- 
LON i i i i i i i ji = j i i lj i j 
96¥ MAX PLOTTED SPEED = 0.20 (M/SEC) sow 


68 


BEST COPY AVAILABLE 


wOPD4 323 02/28/86 


LAYER 1 CURRENTS 


FIGURE 34 
DATE = 349/1966 


os u/s G. OF MEXICO 20132:2: 17.0 
SN T py T T T qe T T T T T en | 
— + . . + + - 
- + ° "> 47 a + + — 
Praa 44% 
Ts hen Bead Dest eazy <vy 33> 
Pneeswereray,)s ftvy + «ev iN 
4byy + <4 <cev¥y 
~ eetil yyy r744245 . + - 
A? >haaneePrPr & a> 
> > 
Sf 7 
_ po ag! “A 


Me 
MY 
ony itt 


4h, 
“, 


S = a Sx at v 
df so a. * ° AX. 
< . 


Seeeadd 


iS ebiel 


~~ 
8 eae 
~ - ~~ 


16N j L l i i les t i i i L 
9eW MAX PLOTTED SPEED = 1.08 (M/SEC) ow 
NORDA 323 02/26/86 
LAYER 2 CURRENTS DATE = 349/1966 
05 u/s G. OF MEXICO 20132:2: 17.0 
JIN T T _s, T T T T a T ae ee T T T 
an + > > + > + > + reer a4 + + + ~ 
a ry veae 444 
< YVa ¥ dene 
«seer = 4479 > > & mad YY Kee e> DDD 
rT << <r ” ° vv Ve caw 
= + + < Yh eennarraaaa enn etatag st Triana ¥ WV be © eediae + + ~ 
*« 4 Va 2>o4 FRP4AA> PT 42 EVE wi «a0 
Arraagvent mA te > > > > F FPP TPT eRe ste 2 4 OO "So © ws 
a<*7 > > 7 P74 Kehoe YO fae a* 
To agga ttt pada tty >» aeittas » 497 F9P, pelee> 2 > benee 4 AZAR A 
4 4 o 
— o)4 48 page POP > EE ANNE MD) Be 4A PEK ee > «6 A + + ~ 
aqe >, eee a > My VOD? & > + o> eke han 4h? > 2 AVY 4G 4M Saal 
4 nea? PtP97 44 TO hc Gdns >> «>> >> a? >» ¥ >e «4 
A@ 14nd J Pewee VE wih de Ceres V4 A Aare sere? 944 4>> 774 Y44 4 4 > ode 
q a anne TA2Ae MAz>D 277d 4 nihydacy <cene © news 944 4A 7 5, > 
v a> > > Pee Ade < < + wa >> saa 
+ 3 bya, de TEA a> > ooaes 40k deer eth el” 7% 448 * ° . 
> Api) nee Taredees, Vin < Far 4 PPT < £ORRA oe ee 4 ght 4 eee a 
th, a Pie ty Bab beeen ALY aqme & Kn Aha 2 44* ¥en > nc 
bah Bf RAVES ond eee enes 25 PeeEy > > Fb ae a) eee ee 
— rt . we > taba 9 47 ane’ t + 
a” Sweet preset PMARR EDT 944 AKET € Ws Ry ‘4, 
> orc 
‘3 Shs exeraha Verte eeeceee «ee cer Pe rig 4 wv 
440° atohetinings rAd Veet $e > 5 qvee ~ “4evwr Y © | 
m feat Fmbhseand *¢ 22 *Asannete cee Y ve ® ee ee ~ 
Ae sss. ave > %s* &« * Y dee? whew vv my 
| See “44 <4 ¢ Pet 4 <$< - avy dae ~ 2. oe! Vvv rn Monee de 
ee ~ ae a2hvV4an Vo» a . rr<* 7 4 - pw tee rr e*eas a? 
m - tee < meer f° re » - abe fede: 7 > “> frweies 
’ Thr Par * 
Ae ahect 44 bee nhs <A A ocioy <v¥»* 2 
a we ae Page Y cen Deve . > ss >> > perere 
i a’ 44 vesqp\ve ee Raw reer nn Vs > A aoaD > =eee" | 
7 7 Phe 
ate ’ eebiv, »YVY VV 
‘..% +* 4<av ‘ Pa 
$80 xy, »<* bee een pie vv A 
. ad «+ 
sites hime vy > © Meee Cee coh v Wy 4 
e + > > <4 <~ * « a 
eet waere ry * . 7944 y4v Ee & dna naz pe eis > > &) 
. AA 7 
rt. «3% Addey vooat ne a y & Hon aaee 
aie hr Ys ves as a’ 
mo reg” >a. v4 in » ¥ #44 . > = 
> are bn DNS 9” . 4 
> ea? See e<*ere Ye Va A 
wot reve fea‘ a er hb 4 CG 
¢3¥ he eee, + PFA van, aA 
b— . sec en 4 Fey ew me ‘a 4 + + + 
Rs inn oe? A ewasal Wea A 
4 24 G5" Debian ey Ve 4 
PZ ov Ase ebbees bh 9 
> mem veeras s + J 
= s Pheee Ce & & VAD 44 + + + 
acces thay | 
<c< smay ? | 
- } 
10M Ld — j 1 i i —— a a 
96¥ MAX PLOTTED SPEED = 0.27 (M/SEC) aow 


BEST COPY AVAILABLE 


wOPO4 323 02/28/86 


69 


FIGURE 35: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/16.0 on model days (a) 2880, 
and (b) 2970. The assigned dates (344/1967 and 069/1968) 
indicate the applied wind forcing, but the experiment was 
not a hindcast and the ccean currents in the Gulf on that 
date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, 
i.e., every 0.2 degrees, and all velocities greater than 50 
em/sec are plotted as 50 cm/sec. 


FIGURE 36: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/16.0 on model days (a) 3060, 
and (b) 3150. The assigned dates’ (1159/1968 and 249/1968) 
indicate the applied wind forcing, but thé experiment Was 
not a hindcast and the ocean currents in the Gulf on that 
date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, 
i.e., every 0.2 degrees, and all velocities greater than 50 
em/sec are plotted as 50 cm/sec. 


FIGURE 37: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/16.0 on model days (a) 3240, - 
and (b) 3330. The assigned dates’ (339/1968 and 0063/1969) 
indicate the applied wind forcing, but the experiment was 
not a hindcast and the ocean currents in the Gul? on that 
date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, 
i.e., every 0.2 degrees, and all velocities greater than 50 
cm/sec are plotted as 50 cm/sec. 


FIGURE 38: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/16.0 on model days (a) 3420, 
and (b) 2510. The assigned dates (153/1969 and 243/1969) 
indicate the applied wind forcing, but the experiment was 
not a hindcast and the ocean currents in the Gulf on that 
date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, 
i.e., every 0.2 degrees, and all velocities greater than 50 
cm/sec are plotted as 50 cm/sec. 


FIGURE 39: Instantaneous view of the upper layer averaged 
velocities from Experiment 201/16.0 on model days (a) 3600, 
and (b) 3690. The assigned dates (333/1969 and 058/1970) 
indicate the applied wind forcing but thé experiment ‘as 
not a hindcast and the ocean currents in the Gulf on that 
date might have been quite different from those shown. 
Vectors are only plotted at every second model grid point, 
i,e., every 0.2 degrees, and all velocities greater than 50 
cm/sec are plotted as 50 cm/sec. 
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FIGURE 40: Location of moorings deployed in the eastern 
Gulf of Mexico for years 1 and 2 of the MMS sponsored Gulf 
of Mexico Physical Oceanography Study, from Waddell 
(1986). Locations V, W, X, Y, 2 in the ocean model 
approximately coorespond to moorings G, A, C, D, F 
respectively. 


PIGURE 41: Kinetic energy spectra at 100 m depth for 
(solid) experiment 201/13.0 at location X, and (dashed) 
mooring C. Variance is shown on a linear scale, but a 
logarithmic scale would be appropriate. 


PIGURE 42: Kinetic energy spectra at 100 m depth for 
(solid) experiment 201/13.0 at location Z, and (dashed) 
mooring F. Variance is shown on a linear scale, but a 
logarithmic scale would be more appropriate. 
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FIGURE 43: Instantaneous view of the layer depth at hour 
0 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topography. The layer depth is positive, but 
very close to zero in the second basin. 


FIGURE 44: Instantaneous view of the layer depth at hour 
6 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topograp.,y. Water is being added to the center 
of the first basin. 


FIGURE 45: Instantaneous view of the layer depth at hour 
12 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topography. Water is being added to the center 
of the first basin. 


FIGURE 46: Instantaneous view of the layer depth at hour 
18 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topography. Water is being added to the center 
of the first basin. 


FIGURE 47: Instantaneous view of the layer depth at hour 
24 from a simulation using a 2-D explicit version of the 

NORDA/JAYCOR ocean model that allows layer interfaces to 

intersect topography. Water is being added to the center 
of the first basin. 


FIGURE 48: Instantaneous view of the layer depth at hour 
30 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topography. Water is being added to the center 
of the first basin. 


PIGURE 49: Instantaneous view of the layer depth at hour 
36 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topography. Water is being added to the center 
of the first basin. 


FIGURE 50: Instantaneous view of the layer depth at hour 
42 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topography. Water is beiig added to the center 
of the first basin. 


FIGURE 51: Instantaneous view of the layer depth at hour 
48 from a simulation using a 2-D explicit version of the 
NORDA/JAYCOR ocean model that allows layer interfaces to 
intersect topography. Water is being added to the center 
of the first basin. 
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FIGURE 44 
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FIGURE 45 
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FIGURE 46 
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APPENDIx A 


APPENDIX A: MODEL ae on \ 
{ 


YEAR 1 REFERENCE PARAMETERS (EXPERIMENT 40): 


upper layer inflow transport = 20 x 10© m3 sec-l 
(= 20 Sv), 


lower layer inflow transport = 10 x 106 m3 sec~l 
(= 10 Sv), 


wind stress = 0, 

horizontal eddy viscosity, A = 300 m2 sec7l, 
grid spacing = 20 by 22 km (0.2 by 0.2 degrees), 
upper layer reference thickness, Hl = 200m, 
lower layer reference thickness, H2 = 3300m, 
minimum depth of bottom topography = 500m, 

beta, df/dy = 2 x 10711 m-1! sec-!l, 

Sos)esse pareneter at the southern boundary, £ = 4.5 x 
10 sec™+, 

gravitational acceleration, g = 9.8 m sec~2, 
reduced gravity, g' = .03(Hl1 + H2)/H2 m sec™2, 
interfacial stress = 0, 

coefficient of quadratic bottom stress = .003, 


time step = 1.5 hours. 
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EXPERIMENT 9: 


26 x 106 m3 sec-l 


upper layer inflow transport 
(= 26 Sv), 


4 x 10© m3 sec-l 


lower layer inflow transport 
(= 4 Sv), 


grid spacing = 25 by 25 kn, 
lower layer reference thickness, H2 = 3400m, 


Coriolis parameter at the southern boundary, £ = 5 x 
10-5 sec-!l, 


coefficient of quadratic bottom stress = .002, 
time step = 1 hour, 


all other parameters as in the Year 1 reference 
experiment. 


EXPERIMENT 31: 


upper layer inflow transport = 0, 
lower layer inflow transport = 0, 


Wind stress from seasonal climatology based on ship 
observations, 


coefficient of quadratic bottom stress = .002, 


all other parameters as in the Year 1 reference 
experiment. 


EXPERIMENT 40: 


the Year 1 reference experiment. 


EXPERIMENT 68: 


Wind stress from 12 hourly Navy Corrected Geostrophic 


Wind set, 


all other parameters as in the Year 1 reference 
experiment. 
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YEAR 2 REFERENCE PARAMETERS (EXPERIMENT 201/13.0): 


20 x 10 m3 sec-l 


upper layer inflow transport 
(= 20 Sv); 


lower layer inflow transport = 10 x 106 m3 sec™l 
(= 10 Sv), 


wind stress = 0, 

horizontal eddy viscosity, A = 100 m2 sec~l, 
grid spacing = 10 by 11 km (0.1 by 0.1 degrees), 
upper layer reference thickness, Hl = 200m 

lower layer reference thickness, H2 = 3450m, 
minimum depth of bottom topography = 500m, 

beta, df/dy = 2 x 1071! m-1 sec-l, 


Coriolis parameter at the southern boundary, f = 4.5 x 
1075 sec -l, 


gravitational acceleration, g = 9.8 m sec™2, 


reduced gravity, g' = .03 m sec~2 
interfacial stress = 0, 
coefficient of quadratic bottom stress = .002, 


time step = 1.0 hours. 


EXPERIMENT 202/11.0 


@® upper layer inflow transport = 0, 


lower layer inflow transport = 0, 

Wind stress is monthly climatology for 6 years and then 
monthly winds for 1967, 1968, 1969, etc. from the Navy 
Corrected Geostrophic Wind data set. 

horizontal eddy viscosity, A = 300 m2 sec~l, 

grid spacing = 20 by 22 km (0.2 by 0.2 degrees), 


time step = 1.5 hours, 


all other parameters as in the Year 2 reference 
experiment. 
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EXPERIMENT 201/13.0 


the Year 2 reference experiment. 


EXPERIMENT 201/16.0 


@ wind stress is zero for 6 years and then monthly winds 


for 1967, 1968, 1969, etc. from the Navy Corrected 
Geostrophic Wind data set. 


interfacial stress is space- and layer- 
thickness-dependent: for an upper layer thickness of 
200 m the frictional e-folding time varies from 3 years 
in the east to 115 days in the west of the Gulf, for a 
layer thickness of 100 m the range is 75 days to 7 
days, and for a layer thickness of 50 mit is 4.5 to 
0.5 days. 


time step = 0.5 hours, 


all other parameters as in the Year 2 reference 
experiment. 


EXPERIMENT 201/17.0 


horizontal eddy viscosity, A = 50 m2 sec~l, 
coefficient of quadratic bottom stress = .003, 


all other parameters as in the Year 2 reference 
experiment. 
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